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The dependence of the g-factors of semiconductor donors on applied electric and magnetic fields 
is of immense importance in spin based quantum computation and in semiconductor spintronics. 
The donor g-factor Stark shift is sensitive to the orientation of the electric and magnetic helds and 
strongly influenced by the band-structure and spin-orbit interactions of the host. Using a multi- 
million atom tight-binding framework the spin-orbit Stark parameters are computed for donors in 
multi-valley semiconductors, silicon and germanium. Comparison with limited experimental data 
shows good agreement for a donor in silicon. Results for gate induced transition from 3D to 2D wave 
function confinement show that the corresponding g-factor shift in Si is experimentally observable. 

PACS numbers: 71.55.Cn, 03.67.Lx, 85.35.Gv, 71.70.Ej 



Understanding the behavior of single donor electron 
bound states under mcsoscopic electric and magnetic 
fields is a fundamental issue critical to current minia- 
turization of semiconductor devices [l[ and to the devel- 
opment of new quantum technologies 

UMM- It is only 

very recently that convergence between experiment and 
theory has occurred for the electric gate control of the 
orbital states [l[ and electron-nuclear hyperfme interac- 
tion [B], Q for a donor in silicon. However, the Stark 
shift of the donor spin-orbit interaction, which is central 
to understanding the precise spin properties in combined 
electric and magnetic fields, is just beginning to be under- 
stood ■ We report the first atomistic treatment of the 
donor spin-orbit interaction in multi-valley semiconduc- 
tors in gated environments, and show non-trivial agree- 
ment with experiment where available. We calculate the 
donor g-factor shift for the transition from 3D Coulomb 
to 2D interface confinement and show that the effect is 
experimentally observable. 

Wave function engineering of donor spins is a basic in- 
gredient of several quantum computing schemes @, S, 0] ; 
and may also help realize novel devices based on spin 
degrees of freedom. In one method, an applied E-field 
deforms the donor wave function and modifies its orbital 
angular momentum, which in turn can modify its spin 
properties through the spin-orbit interaction. This spin- 
orbit Stark effect is manifested by an E-ficld dependence 
of the effective g-factor, and can be probed by ESR ex- 
periments However, as we are dealing with donor lev- 
els in the solid-state, generally with complicated multi- 
valley orbital-spin effects, the physical origins of this phe- 
nomenon is at present not well understood. In semicon- 
ductors with significant spin-orbit interaction, this tech- 
nique can even provide a way to rotate spins by electri- 
cal modulation of the g-tensor, and was demonstrated in 
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FIG. 1: Wave function and g-factor engineering of donor and 
interface confined electrons by EM fields in Si and Ge semi- 
conductors (CB equivalent energy surfaces shown left). 

GaAs quantum dots @ and in GaAs/Al^Gai-^As het- 
erostructures 0, [l(| • While the spin-orbit interaction in 
Si is relatively small, in a quantum computer application 
such effects can lead to qubit errors at the threshold level 
and need to be characterized and understood. 

In this letter we report the first investigation of the 
Stark shift of the donor g-factor in two multi- valley semi- 
conductors with varying degrees of spin-orbit interaction. 
Wc quantify the way in which the E-field removes the 
isotropy of the donor g-tensor components, resulting in 
an anisotropic Zeeman interaction. We also compare our 
g-factor Stark shift for P donors in Si against the only 
available measured value Q, and report corresponding 
parameters for Ge host under different orientations of E 
and B fields, as a guide for future experiments. Finally, 
we investigate g-factors of donors close to an oxide semi- 
conductor interface, and study the g-factor variation as 
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the electron undergoes a symmetry transition from 3D 
Coulomb to 2D interfacial confinement [![. This transi- 
tion is central to proposals for donor-gate-confined inter- 



facial transport and qubits in Si [11 , 

Engineering the magnetic field response in semicon- 
ductors typically involves compound structures such as 
Al^Gai-^As or Si^Gei-^ with a spatially varying mate- 
rial composition. Since the two materials, Al and Ga in 
Al^Gai-^As for example, have different g-factors, the ef- 
fective g-factor of an electronic wave function can be con- 
trolled by pulling the wavefunction from an Al-rich part 
of the device to the Ga-rich part by means of an E-field 
The direct dependence of the g-factor on the field, 
however, has been largely ignored in literature except for 
Refs [l(|, where the g-tensor modulation resonance was 
used in Al^Gai-^As hetero-structures to control spin co- 
herence electrically. In Rcf [lj], 



an all electrical control 
of the spin of a Mn hole in GaAs was investigated. A 
large anisotropic Zeeman splitting has been reported for 
acceptor levels in SiGeSi quantum wells [1 51] and also for 
single quantum states of nanoparticles [16| . Past ESR ex- 
periments 
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18(| have investigated the effect of uniaxial 
strain on donor g-factors in Si and Ge, while a recent 
work demonstrated the gate control of spin-orbit inter- 
action in a GaAs/AlGaAs quantum well 19j |. 

In this work, we employed an atomistic tight-binding 
(TB) theory with a 20 orbital sp 3 d 5 s* basis per atom 
including nearest neighbor and spin-orbit (SO) interac- 
tions. The total Hamiltonian of the host and the donor 
under an applied E-field can be expressed as, 

H = H - -^T^S -pX Wo + t/do„orM + eE ■ f. (1) 

The first term represents the host semiconductor, the 
second term the SO interaction of the host due to the 
crystal potential Vq, the third and the fourth represent 
the donor potential and the applied E-field. The semi- 
empirical TB parameters for Si and Ge 2fJ, 21 1 used here 
have been well established in literature. The SO inter- 
action of the host was represented as a matrix element 
between the p orbitals of the same atom after Chadi , 
and has been shown to cause energy splitting between 
the split-off hole (SH) band and the degenerate manifold 
of the light (LH) and heavy (HH) hole bands. This rep- 
resentation includes both the Rashba and Dresselhaus 
terms inherently, as opposed to the k-p method where 
the two are separately expressed. The donors are repre- 
sented by a Coulomb potential screened by the dielectric 
constant of the host. The potential at the donor site Uq 
was adjusted to obtain the ground state (GS) binding en- 
ergy [23j taking into account the valley-orbit interaction 
in multi- valley semiconductors (24J. The total Hamilto- 
nian was solved by a parallel Block Lanczos algorithm 
to obtain the relevant donor states. A typical simula- 
tion involved about 3 million atoms, and requires about 
5 hours on 40 processors. The Zccman Hamiltonian was 



then evaluated perturbativcly, using the matrix elements 
H Zi j = (y i (r,E)\(L + 2S)-B\V J (r,E)), where i,j rep- 
resent | and J, spins of a donor state. The small B-ficld 
(IT) used throughout this work justified the inclusion of 
linear B-field dependencies only in the Zeeman Hamilto- 
nian. The g-factor was then evaluated using the lowest 
spin states (e), g(E,B) = (e T - e{]/n B \B\. 

This TB model has been previously used to investigate 
the Stark shift of the hyperfine constant for a P donor 
in Si Q in good agreement with ESR measurements [f| 



and momentum space methods 13]. It has also been 
successfully applied to interpret orbital Stark shift mea- 
surements on single As donors in Si FinFETS [l|]. 

The g-factor for a donor ground state in a multi- valley 
semiconductor is influenced by two main factors. Within 
a single valley, the g-factor of electrons moving paral- 
lel to the valley axis (<?||) is different from the g-factor 
for perpendicular motion (<?h), assuming the semiconduc- 
tor has non-spherical energy surfaces. This anisotropy 
may be affected further by external perturbations such as 
strain or E-fields, which may cause higher lyin g co nduc- 
tion bands (CB) to admix with the lowest CB [18[. Sec- 
ondly, the donor ground states in Si and Ge have an equal 
admixture of all the valleys due to the valley-orbit inter- 
action [ii]], resulting in an isotropic effective g-factor. 
Since an E-field removes the equivalency of the valleys, 
the effective g-factor becomes anisotropic depending on 
the contribution of the different valleys to the quantum 
state. The first effect is termed as the single-valley ef- 
fect while the second the valley repopulation effect [l7j |. 
In a tight-binding description, it is not necessary to sin- 
gle out different valley contributions since the full band 
structure is considered in the formalism. Hence both ef- 
fects are captured in the resulting g-factor. 

TABLE I: Comparison of the quadratic g-factor Stark shift 
coefficient for donors in Si and Ge under different orientation 
of electric and magnetic fields. 



Donor 
(valleys) 
[direction] 


Valence band 
splitting 
[eV] 


Binding 
energy 
[meV] 


E-ficld 


B-ficld 


V2 [10" 3 Aim 2 /V 2 ] 


Theory 


Expt [5j 


Si:P 
(6) [100] 


0.044 


-45.6 


[010] 


B|| 


-0.012 


-0.01 
(Si:Sb) 


B± 


0.014 


Gc:P 
(4) [111] 


0.29 


-12.8 


[010] 


B\\ 


-4.8 




B± 


-4.8 


[111] 


B\\ 


143.8 


B ± 


-80.1 



Table [J compares the spin-orbit properties of donors in 
Si and Ge. The SO interaction is stronger in Ge than in 
Si, as shown by the spin-orbit interaction whose strength 
is shown by the energy splitting of the split-off valence 
band from the degenerate light and heavy hole bands at 
the gamma point of its band structure. Si and Ge are 
multi-valley semiconductors with valleys located along 
[100] and [111] crystal axes respectively. 

FiglHshows that the g-factor of donors primarily varies 
quadratically with the E-field. The g-factor shifts are af- 
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FIG. 2: Relative change in the the donor g-factor in Si and 
Ge as a function of E-field strength (in applied direction) and 
for parallel and perpendicular magnetic field orientation, a) 
Si:P under [010] E-filed , b) Ge:P under [010] E-field, c) Ge:P 
under [111] E-field, and d) Ge:P under [110] E-field. The g- 
factor shift is sensitive to the relative orientation of the E and 
B fields with respect to the valley axis. 

fected by the relative angles between the E-field, the val- 
ley axis, and the B-field. Fig [2^i shows the g-factor of 
a P donor in Si subjected to [010] E-ficlds, while b, 
c, and d arc for a P donor in Ge under various orienta- 
tions of the E-field. In Si the [010] E-field (2a) removes 
the equivalency of the six valleys, introducing the vallcy- 
repopulation effect in the donor wavefunction. This re- 
sults in two different parabolas for the g-factor shifts, one 
for B parallel to E and the other for perpendicular ori- 
entation. The [010] directed E-field cannot remove the 
equivalency of the [111] valleys in Ge, and both paral- 
lel and perpendicular B-fields produce the same g-factor 
shifts (03). However, when the field is directed along 
the [111] valley axis as in 2c, we obtain the split g-factor 
parabolas §Bp) similar to Fig 2a. In the absence of an 
E-field, the Zeeman effect of the donor ground state is 
isotropic as shown by the convergence of the two parabo- 
las at E = in both Figs Hi and[2fc. 

The results of Fig [2] are fitted to a quadratic equation 
g(E)/g(Q) — 1 = JfcE 2 , where rj2 is the quadratic Stark 
coefficient. Values of 772 are shown in Table [T] for a few 
different E and B field orientations. For donors close to 
interfaces, there can be some linear Stark effect as well, 
which we discuss briefly in Fig 4. The quadratic coeffi- 
cient of Si:P is of the order of 10~ 5 /um 2 /V 2 , and com- 
pares well in magnitude with the only available experi- 
mental data on Si:Sb. Order of magnitude comparison 
of 772 for Si and Ge shows that the spin-orbit Stark ef- 
fect is stronger in Ge than in Si. The Zeeman anisotropy 
is also stronger for donors in Ge, where the quadratic 
coefficients can differ by an order of magnitude between 
parallel and perpendicular B-ficlds (Table l).The direc- 
tion of the E-field relative to the valley-axes also affects 



the strength of the Zeeman interaction. This is shown 
by comparing the Bu results of Ge:P for [010] and [111] 
directed E-fields. The quadratic Stark coefficients in this 
case differ by two orders of magnitude. Our results show 
very good agreement in magnitude of 772 for Si:P with 
the measured value reported in Ref however, there is 
a discrepancy regarding the relative sign of the g-factor 
shift in the quadratic regime which we are currently un- 
able to account for. Within the TB framework calcula- 
tions of the g-factor Stark shift for the single valley GaAs 
case were also carried out and showed good agreement in 
magnitude in comparison with the recent k.p results Q. 

A simple multi-valley picture provides some intuitive 
explanations of the Stark shifted g-factor based on the 
valley repopulation effect. If \a x \ 2 represents the contri- 
bution of the +x valley in Si to the donor GS, and g, x 
the diagonal g-tensor corresponding to this valley with 
the x component given by g\\ while the y and z compo- 
nents given by g±_ , then the effective g-tensor of the donor 
GS is given by g(E) = T, i= ± x ,± y ,± z \ a i?9i- Assuming 
dj = a_; and a x = a z for a [010] E-field, we obtain the 
effective g-tensor components g x , g y and g z as, 



9x 
9y 



9z = 

41 GU 



2(\a y \ 2 + \a x \ 2 )g ± + 2\a x \ 2 gil 



9\\ 



(2) 
(3) 



These equations show that the parallel component of 
the g-factor has a different response to the electric field 
as compared to the perpendicular component, verifying 
the split g-factor curves of Fig [3^ and [5U. At E = 0, 
each di = l/v 7 ^, and eqs (2) and (3) reduce to g x = g y = 
g z = !<jLi_ + |<7|| = g , showing an isotropic effective g- 
factor. At ionizing E-fields, only the valleys parallel to 
the field contribute to the state. Setting a y = \/y/2 and 
o> x = o-z = in (2) and (3), g x = g± and g y =511, which 
helps to probe the single valley g-factors, as shown later 
in Fig|4] Similar expressions can be derived for a donor in 
Ge taking into account that the Ge valleys are in [111]. 
For a more quantitative approach, however, one needs 
to know also the g-factor variation within a single val- 
ley, the precise nature of the wavefunction distortion by 
the E-field, and the effect of B-ficlds. The TB approach 
provides a generalized framework to include all these. 

In Fig [3J we vary the angle 6 between the E and B- 
fields from to 90 degrees for a) Si:P under [010] E-field, 
and b) GeP under [111] E-field. The relative change in 
g-factor shows a linear dependence on sin 8, consistent 
with Ref [TH. The sensitivity of this variation increases 
at higher E-fields as shown in Fig [5^ and [SJs. The flat 
E = line indicates that the Zeeman effect is isotropic at 
zero field. Assuming an y-directed E-field in Si such that 
the effective g-tensor diagonal components are given by 
g x {E) = g z {E) ^ g y {E), the linear dependence of g(E) 
on sin 2 9 is shown by, 
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FIG. 3: Anisotropic Zeeman effect with E-field. 9 represents 
the angle between the E and the IT B field, a) Si:P under 
[010] E-field with B-field varying from [010] to [001]. b) Ge:P 
under [111] E-field, with B-field varying from [111] to [1-10]. 

Eq 2] can be derived by expanding g — (g y (E) 2 cos 2 + 
g x (E) 2 sin 2 0) 1 / 2 up to linear terms in (g x (E) 2 / g y (E) 2 — 
1) sin 2 0. 
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function becomes asymmetric due to sharp truncation by 
the surface. Fig [4] shows at small donor depths, the g\\ 
and g± also exhibit a linear field dependence. In this 
regime, the linear Stark coefficient exceeds the quadratic 
coefficient. A similar effect was observed for the hyper- 
fine Stark effect @. 

In conclusion, we have for the first time applied atom- 
istic techniques to understand and predict the E-field re- 
sponse of donor g-factors in multi-valley (Si, Ge) semi- 
conductors with different degrees of spin-orbit interac- 
tion. The E-field induces a Zeeman anisotropy that varies 
with the relative angle between the E and the B fields. 
The strength of the Stark shift is also dependent on the 
direction of the E-field relative to the valley axis. The 
computed Stark shift coefficient of Si:P compares well 
in magnitude with the only available measured value for 
donors in Si. The donor g-factor Stark shift was also 
computed for the 3D to 2D confinement transition, sug- 
gesting that the effect is accessible to experiments. 
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FIG. 4: Interface effects on the donor g-factor for Si:P. g- 
component a) parallel, b) perpendicular to the field and the 
interface at various donor depths. 

The confinement transition from 3D Coulomb to an 
interface 2D system has recently been observed (in con- 
junction with the theoretical approach used here) 
and is related to a new control scheme based on dopants 
close to the Si-Si02 interface llj. The donor electron 



can be adiabatically pulled to the interface by gate volt- 
ages [lU dll [2^], and controlled by surface gates. We 
computed the g-factor of a system undergoing this con- 
finement transition. Fig 0] shows the components of the 
g-factor parallel (|3Ji) perpendicular (0J)) to the E-field for 
various donor depths. As the E-field increases, the two Si 
conduction band valleys in the direction of the field are 
lowered in energy relative to the four valleys perpendic- 
ular to the field axis. The interface state realized at ion- 
izing E-ficlds has contribution from these two uni-axial 
valleys, and hence their gu and g± approach those of the 
two valleys. Our simulations indicate gu — g± ps 8 x 10 -3 , 
which compares well in order of magnitude with the mea- 
surements of Ref [l8| . T his g -factor anisotropy has also 
been reported in 2DEGS [27J ■ The transition to the sin- 



gle valley g components is abrupt if the donor is far away 
from the interface (>10 nm ), and g radual if the donors 
are closer to the interface [13J, [25|, |26( . Proximity to inter- 



faces is also marked by linear Stark effect since the wave 
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